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The anhamionic vibrational dynamics in nearly stoichiometric La2Cu04+5 is studied by means 
of anelastic and ^^^La NQR relaxation. In the absorption component of the elastic susceptibility as 



o 

fj ■ well as in the nuclear relaxation rate a peak is detected as a function of temperature, and a relaxation 

I ' time T = 1.7 x 10^^^ exp[(2800 K)/T] s is derived. The relaxation processes are attributed to tilt 

^\, motion of the CuOe octahedra in doublewell potentials, whose cooperative character increases the 

r^ ' effective energy barrier to the observed value. The analysis of the relaxation mechanisms has been 

C/3 ' carried out by reducing the dynamics of the interacting octahedra to a one-dimensional equation 

of motion. The soliton-like solutions correspond to parallel walls separating domains of different 
tilt patterns and give rise to pseudo-diffusive modes which appear as a central component in the 
spectral density of the motion of the octahedra. The tilt waves may be considered to correspond to 
the dynamical lattice stripes observed in La-based and Bi-based high-Tc superconductors. 



I. INTRODUCTION 

There is growing interest in the local deviations from 
the average structure recently observed in cuprate su- 
perconductors. These local inhomogeneities are of- 
ten identified with one-dimensional conducting domains 
of increased carrier density and antiferromagnetic car- 
rier depleted domains.EI In some theoretical approaches 
such stripes are believed to play a role in the mech- 
anism underlying high-Tc superconductivity. Recently, 
the anelastic relaxation spectrum of undoped La2Cu04 
has been found to exhibit two intense relaxation pro- 
cesses, attributed to intrinsic relaxational dynamics of 
unstable lattice modes, namely quasi-rigid tilts of the O 
octahedra.El The relaxation process below 30 K is related 
to an instability toward the LTT structure, while the 
high-temperature process is thermally activated with an 
activation energy around 2800 K. In the present report 
combined anelastic and ^^^La NQR relaxation measure- 
ments of La2Cu04+5 are presented and interpreted in 
terms of the collective dynamics of tilting CuOg octa- 
hedra in double- well potentials. The tilt waves can be 
considered the equivalent in the undoped system of the 
lattice stripes and structural inhomogeneities observed in 
high-Tc superconducting oxides. 



II. EXPERIMENTAL AND RESULTS 

The sample was a sintered bar 45 x 4.5 x 0.5 mm'^ of 
La2Cu04+5; the excess O was extracted by outgassing 
at 750 °C down to an O2 partial pressure of the order of 
10~* torr. Such treatments can introduce few at% of O 
vacancies in the Cu02 planesjj which are reversibly filled 
after reoxidizing at high temperature. These vacancies 
do not affect the anelastic or the NQR relaxation, being 
static at the temperatures of interest (below 500 K). The. 
O stoichiometry in the as-prepared state is estimatedtl 
6 ^ 0.004. Details on the sample preparation and char- 
acterization can be found in Ref. |[ 

The elastic energy loss coefficient Q^^, measured 
by exciting the flexural vibrations, is given by 
S" (w) /S' (w) ~ S" (w) /S (00), where S {u) is the com- 
plex dynamic compliance, with unrelaxed value S'(oo). 
The contribution to the imaginary susceptibility from a 
relaxation process with characteristic time r isci 



-^oc [x{SXf /T] LOT/[l + {u;Tf 



(1) 



where x is the atomic fraction of relaxing entities, each 
producing a change ^A of the strain e when its state 
changes; uj/2tt is the measuring frequency. The Q~^ (T) 
curve has a maximum at the temperature where ujt (T) — 
1, yielding the relaxation rate at the peak temperature. 
The main processes appearing in the anelastic spectrum 
of La2Cu04+5 above 100 K are due to the diffusive jumps 
of the interstitial O atoms. The peak attributed to the 
tilts of the O octahedra appears only after reducing the 



content of excess O down to very low levels {6 < 0.001). 
In fact the interstitial O atoms block ike relaxational 
motion of several surrounding octahedra.a The anelastic 
peak in the nearly stoichiometric state is shown in Fig. 1 
(left hand scale). The two sets of data refer to the exci- 
tation of two flexural modes (460 Hz and 5.9 kHz). The 
peak due to O diffusion was 30 times smaller than in the 
as prepared state, indicating a concentration of residual 
oxygen S ^ 1 x 10^**; such a small contribution has been 
subtracted in Fig. 1. 

The ^^^La NQR relaxation rate Wq in the same sam- 
ple is also shown in Fig. 1 (right hand scale), again 
after subtraction of the contribution from the resid- 
ual interstitial O which gives a monotonous increase of 
Wq for T > 400 K.eI In Lii2-xST^Cu04+s the ^^^Ln 
NQR relaxation, can be due to magnetic or quadrupo- 
lar mechanisms Jj'B The magnetic relaxation is due to the 
fluctuation of the field at the La site, arising from the 
time dependence due to the Cu"'""'" spin operators. The 
quadrupolar mechanism is due to the fiuctuation of the 
transverse components of the electric field gradient at the 
La site, mainly caused by the vibrations of the surround- 
ing O ions, or by a diffusing interstitial O atom passing 
close to the La nucleus. The two tyj>es of relaxation 
are described by different recovery laws,Q so allowing the 
identification of the relaxation mechanism by irradiating 
different NQR lines. The coincidence in Fig. 2 of the 
recovery plots of the NQR echo signals at 2i'q and Svq 
(i/Q = 3el4^ Q/[2I{2I - 1)] « 6.2 MHz, eQ is the ^^SLa 
quadrupole moment) for T > 77 K indicate that above 
this temperature the mechanism is totally quadrupolar, 
i.e. due to the diffusion of interstitial O (subtracted in 
Fig. 1) and to the motion of the apical O atoms. The re- 
laxation rate t~^, determined from the condition ujt = 1 
at the maxima of the Q~^ (T) and Wq (T) curves, is 
plotted in Fig. 3 in logarithmic scale against the recip- 
rocal of temperature. The three points can be closely 
fitted by a straight line t~^ = Tq^ exp{—Ecs/kBT), 
To = 1.9 X 10^12 s ajjd E^s/kB = 2800 K, clearly indicat- 
ing that the same process is observed by both techniques. 



III. DISCUSSION 

The relaxation character of the responses imply that 
the apical O atoms in the LaO planes perform pseudo- 
diffusive motion in a strongly anharmonic potential. 
Such a potential should be associated with the softer lat- 
tice modes of the rigid-unit type, which are low frequency 
modes involving relative rotations of rigid lattice units, 
with little or no distortion of the units (which would in- 
crease the restoring force and frequency of the mode). 
One can assume that the displacement of an apical O 
atom represents the rotation of the rigid octahedron, and 
both the anelastic and NQR relaxation are due to the 
switching of the octahedra between different tilts in mul- 
tiwell potentials. Such potentials have been calculated 



for La2Cu04la and optiaially doped La2_a:Baa;Cu04+5t£l 
and La2-2:Sr2:Cu04+a.lld Minima were found for tilts cor- 
responding to displacements of the apical O atoms by 
S '^ ±0.2 A in the [110] directions (parallel to the in- 
plane CuO bonds), corresponding to the low-temperature 
orthorhombic (LTO) structure, separated by an energy 
barrier of S = 35-53 meV (400-600 K). Additional min- 
ima are found in correspondence to the low-temperature 
tetragonaL.(LTT) structure with tilts about the [100] 
directions,^ but the LTT domains can be observed by 
diffraction experiments only when doping around g by 
Ba or combined Sr and Nd substitution.Ej In the follow- 
ing both the anelastic and NQR measurements are inter- 
preted in terms of cooperative motion of the octahedra 
in such multiwell potentials. 

Let us first consider the ^■^^La NQR relaxation due 
to the atomic displacements s (t) of the apical O atoms 
surrounding the La nucleus. In order to simplify the in- 
terpretation, we will assume uncorrelated motion of such 
atonis. and approximate the relaxation rate Wq in the 
fornxJ 
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(s (t) s (0)} e^'-^'dt = (2) 



where a is a constant of the order of 5/8/ ~ 0.178, 
(s(i)s(O)} is the autocorrelation function of the motion of 
the apical oxygen and J (uj) the corresponding spectral 
density. If s (t) can switch between two values ±S corre- 
sponding to two tilt orientations with probability 1/ (2r), 
then the correlation function is (s(i)s(O)) = (5^e~l*l/'^ and 
the NQR rate becomes 



WQ=AH^2T/[l + {L0Tf]. 



(3) 



where t is the correlation time. 

The imaginary part of the elastic compliance S" (uj) 
can be expressed in terms of the spectral density of t 
autocorrelation function of the macroscopic strain e (t)l 

S" {uj) ^ {ujV / 2kBT) f {e{t)e{0))e-''^*'dt, (4) 



V being the sample volume. In the assumption that e (t) 
is directly related to the displacements s {t) of the O api- 
cal atoms, Eqs. (||) and (^) differ from Eqs. (||) and (||) 
only by a factor co/T and by the constants expressing the 
dependence of strain and NQR frequency on the atomic 
displacements. 

Now it must be explained why the apparent activa- 
tion energy Ecs for the switching rate r~^ is higher 
thajL_the theoretically estimated local potential barrier 
i^.B^Hl This is related to the cooperative character of 
the motion of the octahedra, which has been widely dis- 
cussed by Markiewicz.Ej He considered an LTT ground 
state and observed that, neglecting intraplane interac- 
tions, the octahedra along directions perpendicular to the 



tilt axes are strongly correlated, since they share the O 
atoms which are displaced from the CuO plane, while 
adjacent rows are weakly correlated, since the shared O 
atoms remain in the CuO plane. The system is thus 
reduced to a one-dimensional array of rows of octahe- 
dra, where all the octahedra belonging to the same row 
are tilted according to the same pattern. The resulting 
equations of motion are non-linear and admit solitonic so- 
lutions, which correspond to one-dimensional propagate, 
ing walls between domains with different tilt patterns. E£l 
This analysis allows one to apply the one-dimensional 
models of non-linear lattice dynamics which have Jaeen 
used to describe structural phase-transformations, Ell as 
well as thcj-cprrelated dynamics of off-centre atoms in 
perovskitesJl3 One can consider a potential of the form 
V = ^- —asf + bsj + cSiSi-^i, where each atom moves in 
two minima separated by an energy barrier E — (a^/46) 
with a bilinear coupling to the neighbors. The coupling 
constant c takes into account a cluster average over con- 
figurations. The resulting equation of motion has soli- 
tonic solutions, similar to those found ior the octahedra 
in La2Cu04,Ej and it can be provedEj that the spec- 
tral density J (w) contains both a resonant peak at the 
frequency of vibration in each well, and a central peak 
of pseudo-diffusional character with characteristi c fre- 
quency r~^ = Tq exp (Ees/T) with EeS « 1.75£^A/2S/a, 
To = - \/c/a, where v is the average velocity of propa- 
gation of the soliton-like excitation through the atoms 
spaced by d. These solitonic solutions provide the ob- 
served relaxational contribution to the spectral density 
and susceptibility. The effective potential barrier E^ff is 
increased with respect to the local barrier E by the in- 
teraction between the octahedra. 

In fitting the experimental data with the above model, 
a distribution of the values of c has been introduced, 
corresponding to the distribution in size and shape of 
the regions where the octahedra clusters build up the 
cooperative dynamics. Such regions are limited by the 
excess O atoms, which block the relaxational dynamics. 
The interaction parameter c was distributed according to 
a gaussian, resulting in a distribution of effective energy 
barriers Ecs- 

A feature of the anelastic peak which cannot directly 
be accounted for by the above formulas is the increase 
of its intensity at higher temperature, instead of a de- 
crease as expected from Eq. (|l|). Such a tempera- 
ture dependence is observed in the case of relaxation 
among states which differ in energy by AE > k^T. 
In this case the relaxation strength must be multiplied 
by a factor containing, the product of their equilibrium 
occupation numbers,tJ 4rjin2 = [cosh(Ai?/2fcBr)]~ . 
This causes a decrease of the relaxation strength below 
T ^ AE/kB- The relaxation rate must be multiplied by 
cosh{AE/2kBT), but this correction does not affect siz- 
ably the relaxation curves. The expression o£the NQR 
rate should be modified in the same manner .tl It should 
be remarked that these corrections are strictly valid for 



relaxation between levels without cooperative effects, and 
their extension to the above model of cluster dynamics is 
not obvious. 

The resulting fit is shown in Fig. f as continuous lines. 
The mean values tq = 1.7 x 10^^^ s and E^s = 2800 K 
are the same as deduced from the condition lot = 1 at 
the maxima (Fig. 2). The effective baEriei: EcS is com- 
patible with the theoretically estimatedoHl local barrier 
of "500 K with a mean coupling constant cq — 5a. The 
width of the distribution of c is cr = O.IScq for the NQR 
and 0.25co for the anelastic data (a temperature depen- 
dent width may result from the ordering of interstitial 
O); the resulting distribution in E^s is gaussian with a 
width of - 300 K. The asymmetry energy A£' = 280 K 
is 10 times smaller than E^s and therefore it does not 
affect the overall picture. 

The fit value for the strength of the NQR relaxation 
mechanism is A ■ 5 = lA kHz; A can be estimated in 
the hypothesis of independent displacements of the O 
atoms, assuming dVzz/ds ~ iVzz/d {d ~ 2.7 A is the 
La-0 distance), and for S ^ 0.2 A it yields a relaxation 
strength A • ^ ~ (S^q/d) 6 ~ 8.6 MHz, about a factor 
100 larger than the experimental one. A reduction of 
the strength of the relaxation mechanism may be due to 
the correlation between the displacements of the five api- 
cal O atoms surrounding the La nucleus, giving rise to 
a form factor which multiplies the spectral density. This 
correction, however, can hardly produce a reduction by 
100 times, which is likely to reflect the fraction of octa- 
hedra which are instantaneously changing the tilt orien- 
tation. In a simple picture one can consider independent 
solitonic domain walls (DW) , corresponding to tilt waves 
which propagate through the lattice with velocity v; each 
La nucleus contributes to the NQR relaxation only when 
run over by the tilt wave, so that the relaxation strength 
should be reduced by the fraction of octahedra within 
the domain walls. 

Along the same line, one could attribute the ob- 
served increase of the anelastic relaxation intensity with 
temperature to a temperature dependent volume frac- 
tion of DW, instead of introducing the asymmetry en- 
ergy AE. An increase of DW density, roughly as 
[Tt — T) , is indeed-observed on approaching marten- 
sitic transformations,E2lbut is not sufficient to explain the 
increase of the anelastic relaxation strength by 30% from 
147 to 170 K (Fig. 1, taking into account the 1/T factor 
in Eq. (|l|)). In the present case, the DW exist below 
the transformation from the high temperature tetrago- 
nal structure to the LTO structure at Tt ~ 530 K and 
their density would increase of only 6% on increasing the 
temperature from 147 to 170 K. 

From the close coincidence of the relaxation rate de- 
scribing the anelastic and NQR results, it is deduced that 
the propagation of soliton-like walls is the origin of both 
types of relaxation, and that there is a direct relation- 
ship between the atomic displacement s and macroscopic 
strain e whose spectral densities are probed by the two 
experiments. Such a direct relationship can be assumed 



for the movement of a wall separating two LTO domains 
with opposite signs of the in-plane shear strain e^x — £yy 
(in LTO notation, with a and h directions at 45° with the 
CuO bonds) . The propagation of such DW changes the 
domain sizes and therefore modulates e. 

The spatial arrangement of the domains of tilted oc- 
tahedra cannot be deduced by our data, but the above 
analysis corresponds to an array of propagating paral- 
lel fronts of tili—waves. Similarly, the picture proposed 
by MarkiewiczO for La2_a;Sr3,Cu04+5 is a LTT ground 
state with lowest energy excitations consisting of par- 
allel propagating LTO walls. A closely spaced struc- 
ture of alternating LTT domains and propagating LTO 
walls would form a so-called dynamic Jahn- Teller phase 
with average orthorhombic structure, and conventional 
diffraction measurements would see the average structure 
(the intermediate structure between two LTO domains is 
LTT and viceversa). Such a microstructure closely cor- 
responds to the lattice stripes proposed by Bianconi et 
al. on the basis of EXAFS-Efisults in the same doped ma^. 
terial Lai.85Sro.i5Cu04+5Lj and in Bi2Sr2CaCu208+ir-t 
Stripes have also been observed by neutron diffractionEl 
in Lai.g_a;Ndo.4Sr3;Cu04 with x = 0.12, where they are 
pinned by the formation of a stable LTT structure. Also 
our ^^^La spectra exhibit asymmetric broadening of the 
lines split, by the local magnetic field observed by other 
authors ,E3 supporting the picture of channels where the 
charges doped by the residual excess O are mobile and 
depress the antiferromagnetic order. 



IV. CONCLUSION 

In summary, thermally activated relaxation processes 
due to strongly anharmonic lattice modes have been ev- 
idenced in nearly stoichiometric La2Cu04+5 from the 
elastic susceptibility and ^^^La NQR relaxation rate ver- 
sus temperature, as a low-frequency tail of the spectral 
density of the atomic motions. The pseudo-diffusive mo- 
tions have been identified as the tilt modes of the O octa- 
hedra in a multi-well potential. In the interpretation, we 
have used the fact that the dynamics of the octahedra in 
the ah plane can be cast into a non-linear one-dimensional 
equation, which admits solitonic solutions corresponding 
to propagating walls between domains with different tilt 
patterns. Thus it has been possible to rely on the analysis 
of similar one-dimensional systems, for which the spec- 
tral density contains a pseudo-diffusive central peak of 



idth 



Tg exp {—Ecs/kBT) related to the soliton- 



like lattice excitations. The measured effective barrier 
Ecs/k-Q — 2800 K is about 5 times higher than the bar- 
rier of the theoretically estimated local potential, due to 
the cooperative character of the motion. The alternating 
LTO domains and LTT walls should correspond to the 
lattice stripes observed in superconducting oxides. 
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Fig. 3. Relaxation rate t~^ (T) corresponding to the 
condition lot = 1 at the maximum of the anelastic and 
NQR relaxation curves (Fig. 1). The slope of the straight 
line provides the effective activation energy Ecs, which 
coincides vifith the mean activation energy of the best fit 
in Fig. 1. 



V. FIGURES CAPTIONS 



Fig. 1. ^•^^La NQR relaxation rates Wq (right hand 
scale) and elastic energy loss Q~^ (left hand scale) for a 
sample with S ^ 10^'', after subtraction of the contribu- 
tion from the O diffusion. Solid lines: best fits according 
to the model described in the text. 

Fig. 2. Recovery plots of the amplitude of the free 
induction signal M after saturation of the 3vq NQR line 
(closed symbols and solid lines) at 77 and 297 K. The 
open symbols are the recoveries corresponding to the irra- 
diation of the 2i^Q NQR line. The dashed and dotted lines 
are the curves expected in the case of pure quadrupolar 
and pure magnetic relaxation mechanisms, respectively. 
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